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Within the non-relativistic potential scattering theory, we derive a generahzed version of the 
Liischer formula, which includes three-particle inelastic channels. Faddeev equations in a finite volume 
are discussed in detail. It is proved that, even in the presence of the three-particle intermediate states, 

CLc the discrete spectrum in a finite box is determined by the infinite-volume elements of the scattering 

2h '. S'-matrix up to corrections, exponentially suppressed at large volumes. 

(N . 

^ ; Pacs: 12.38. Gc, ll.lO.St, 11.80. Jy 

Keywords: Resonances in lattice QCD, field theory in a finite volume, Faddeev equations 

(N 

cn 

o 

(N 



% 



1 Introduction 

The nature of the Roper resonance A^(1440) has been an open question for decades. It is quite 
unnatural - at least, from the point of view of the quark models [1, 2] - that its mass turns out 
to be lower than that of the negative-parity ground state A^(1535). Different phenomenological 
approaches have been employed so far to explain such a level ordering. For example, A^(1440) 
was assumed to be a hybrid state with excited gluon field configurations [3,4], a breathing 
mode of the ground state [5], and a five-quark (meson-baryon) state [6]. The properties of the 
Roper resonance have been studied within the Skyrme model [7] and the bag model [8]. 

Up to now, numerous calculations of the excited spectrum of the nucleon [9-24] (see e.g., 
ref. [25] for a recent state-of-art review) have not resolved the issue. In particular, a reverse 
level ordering (the same as in the constituent quark model) between A^(1440) and A^(1535) has 
been reported in some lattice calculations [14] , albeit the situation there is far from being clear. 
It should be also pointed out that a large chiral curvature for the Roper resonance mass is not 
compatible with the findings of ref. [26], where the dependence of the Roper mass on the pion 
mass was investigated within the framework of Chiral Perturbation Theory. According to these 
findings, the level crossing between A^(1440) and A^(1535) states, emerging as a result of chiral 
extrapolation to the small quark masses, does not constitute to a plausible scenario. 

There is, however, one issue that has not been addressed in all above investigations so far. 
The A^(1440) is a resonance and not a stable state which would correspond to an isolated energy 
level in lattice simulations. This, in particular, means that the energy levels measured on the 
lattice will be volume-dependent and the true resonance pole position should be extracted from 
the volume- dependent spectrum. 

The case of the elastic low-lying resonances on the lattice has been investigated in detail. 
Namely, the Liischer formula [27] enables one to uniquely relate the discrete energy levels in 
a finite box to the elastic scattering phase shift in the infinite volume, measured at the same 
energy. This eventually opens the way for the extraction of the parameters of the elastic 
resonances - their masses and widths - in the lattice QCD (for illustration see, e.g., refs. [28, 
29]). 

The case of the inelastic resonances, however, is more complicated. Albeit the Liischer 
approach can be straightforwardly generalized to the case of the coupled two-particle chan- 
nels [30-32], the disentanglement of physical observables becomes a more delicate affair, since 
there is more than one observable at a single energy. In order to circumvent this problem, in 
refs. [32,33] the use of the twisted boundary conditions or asymmetric boxes was advocated 
for the TTvr — KK coupled-channel system. Moreover, in ref. [33] it has been argued that, using 
unitarized ChPT (UChPT) in a finite volume, it is possible first to directly fit the parameters of 
the chiral potential to the energy spectrum measured on the lattice and eventually to determine 
the physical observables (the phase shifts, resonance parameters, etc) from the solution of the 
scattering equations. The feasibility of such a procedure was demonstrated in the example with 
synthetic data [33]. For the further applications of this method, see refs. [34,35]. 

Inelastic resonances, which have a significant decay rate into the three and more particle 
final states, have received much less attention in the literature so far. The Roper resonance 
iV(1440), which decays into the inelastic three-particle channels with an approximately 40% 



probability, is an example of such a system. A priori, one may expect significant finite-volume 
effects in this decay, which can not be evaluated by using the standard Liischer approach. For 
this reason, at the moment it is not clear, whether the reverse level ordering between iV(1440) 
and A^(1535) in lattice simulations, mentioned above, can not be at least partially attributed to 
these finite-volume effects. Putting it differently, one may ask, whether the effect still persists 
for the true resonance positions in the infinite volume. We therefore conclude that it is highly 
desirable to construct a framework that will allow one to systematically calculate the finite- 
volume effects, coming from the tree-body final states. 

Formulating a counterpart of the Liischer approach in a three-body case represents a major 
challenge. For this reason, at the first step, we want to simplify the problem as much as possible. 
Namely, we consider a non-relativist ic quantum-mechanical model with coupled two-particle 
and three-particle channels. Multi-particle channels (4 and more particles) are neglected from 
the beginning. Using these approximations, in particular, the technical complications, related 
to the necessity of Lorentz-boosting the two-particle sub-systems in the three-particle state 
to their respective center-of-mass (CM) frames can be avoided. Moreover, using the potential 
model instead of the effective field theory (EFT) framework, we avoid the discussion of the 
proper counting rules for the multi-particle intermediate states, which are generated by loops. 
The core of the problem, which consists in the study of the finite- volume effects coming from the 
three-particle intermediate states, remains however unaffected by the above approximations. 
At the next step, we plan to carry out a full-fledged investigation of the problem. 

The central issue addressed in our study, can be formulated as follows. In case of the elastic 
two-body scattering, the Liischer formula relates the discrete spectrum in a finite box to the 
scattering S'-matrix element in the infinite volume, up to the corrections that are exponentially 
suppressed in the box size. Putting it differently, this is a relation between the observables in 
a finite and in the infinite volumes: the details of the potential do not matter if the box size 
is much larger than the typical radius of interaction. It is intuitively clear that the statement 
must remain valid even if there is a contribution from the three-particle intermediate state - as 
far as the size of the box remains much larger than the scales characterizing the interactions. 
In our paper, we prove this statement for the model described in the previous paragraph. It is 
expected that this proof can be extended beyond this particular model. 

The layout of the present paper is as follows. In section 2 we discuss the quantum- mechanical 
model which describes the scattering in the coupled two- and three- particle channels. Con- 
sidering the same model in a finite volume, we derive an equation which predicts the volume- 
dependent energy levels. Further, in section 3 we re-derive the Liischer formula for the two-body 
scattering applying a new method, which can be used with minor modifications in the 3-body 
case as well. This is done in section 4, in which we derive the three-body analog of the Liischer 
formula. Section 5 contains our conclusions. 



2 Three-body problem in a finite volume 

2.1 The model in the infinite volume 

In this section, we describe a model which will be later used to study the finite-volume effects 
in the three-body sector. In order to make the presentation self-contained, below we explicitly 
display the standard formulas in case of the infinite volume. We consider a non-relativistic 
quantum-mechanical system of three spinless non-identical particles with the masses rria, a = 
1,2,3. The Hamiltonian of the model is given by a sum of a free Hamiltonian Hq, pair 
interaction Hamiltonian H2^2 and the Hamiltonian H2^3 that describes the transition from 
two- to three-particle state 1-1-2— )-l + 2-|-3. In order to simplify the following expressions, 
we do not include the three-particle force into the Hamiltonian. As shown later, the induced 
three-particle force will anyway emerge from the two-particle interactions, so one could have 
indeed added it from the beginning without much effort. The explicit expression of the full 
Hamiltonian in terms of the creation/annihilation operators is given by 

H = Ho -|- H2_^2 + H2^3 = Hq -|- H| , 

3 



Ho = ^^(m„ + -— ja„(k)aj,(k), 

^=1 I- V "/ 



a=l k 
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^2^2 = J2 Yl (27r)=^5=^(ki+k2-k3-k4)K;3(ki,k2;k3,k4) 

(Q/3)=(12),(23),(13)ki-k4 

X ac,(ki)a^(k2)4(k3)aJ3(k4), 

H2^3 = ^ (27r)3<5=^(ki + k2-k3-k4-k5) 
ki-ks 

r(ki, ks; k3, k4, kg) ai(ki)a2(k2)al(k3)4(k4)4(k5) + h.c. i , (1) 

where aQ(k) and ajj(k) denote the annihilation and creation operators for the particle a, 
respectively. Further, V12 = V3, V23 = Vi, V31 = V2 denote the pair potentials, so that 
H2^2 = V1+V2+V3 = J2a=i ^a- Finally, the vertex T describes the transition 1+2 — t- 1+2+3. 
Next, we introduce the free and full resolvents, which are defined by 

Go(^) = „ ^ .„ , G(^)= ^ (2) 

Hq — z — zO H — z — zO 

The scattering matrix, which is expressed through the resolvents as 

G(^) = Go(^) + Go(^)T(^)Go(^), (3) 



obeys the Lippmann-Schwinger (LS) equation^ 

T(z) = (-H,) + (-H,)Go(^)T(^). (4) 

In this paper, we are primarily aimed at the modehng of the vrA^ scattering in the Roper 
resonance region - that is, in the presence of an open three-particle channel. To this end, we 
consider the scattering amplitude for the process 1 + 2—7-1 + 2. Introducing the projection 
operators P and Q = I — P (where I stands for the unit operator), which project onto the 
two-particle state |12) and the three-particle state |123), respectively, we obtain the equation 
for the effective two-body operator Tp{z) = PT(z)P 

Tp{z) = W{z) + W{z)Gp{z)Tp{z) , W{z) = (-Hs^s) + H2^3G3(^)(H2^3)^ , (5) 

where the three-particle Green's function obeys the equation 

G3(^) = Gq + Gq(-H2^2)G3(^) , (6) 

and the following notations are used 

Gp = PGo(2)P, Gq = QGo(^)Q. (7) 

Note that the effective one-channel potential W{z) becomes non-Hermitian above the three- 
particle threshold z > mi + m2 + m^. 

Next, we define the three-particle scattering amplitude R(^) through 

G3(^) = Gq + GqR(z)Gq . (8) 

The quantity R{z) can be expressed as 

3 

R(^)= 5^ M,^(^), (9) 

a,/3=l 

where Mq,/3(z) obeys Faddeev equations (see, e.g., [36]) 

3 

M^f,{z) = 6^^T^{z) + T^{z)GQiz) ^(1 - 6^^)M^^{z) , (10) 

7=1 

and Ta{z) denote the scattering amplitudes in the channel a 

T^{z) = (-VJ + (-VJGq(z)T„(z) . (11) 

In order to write down these equations explicitly in momentum space, it is useful to work in 
Jacobi basis 

P = ki + k2 + k3, Pa = K, ^^^ m,kf,-mpk, ^ (a/37) = (123), (231), (312). (12) 

m^ + m^ 



^Note that the sign convention in the LS equation below coincides with the one in the field theory and is 
opposite to the one usually adopted in the potential scattering theory. 



Separating the center-of-mass (CM) motion from the matrix elements 
{k[k',k',\M^f,{z)\k,k2ks) = {27ifS%k[ + k'2 + k:, - ki - k2 - k3)(p:,q:,|m,/3(z)|p^q/3) 

(k;k^k'3|GQ(2)|kik2k3) = (27r)353(k; + k^ + k;,-ki-k2-k 



X 



3; 

2 2 ' 

2Ma 2/i„ 

(k;k^kyT,(2)|kik2k3) = (27r)353(k; + k^ + k'3 - ki - k2 - k3)(27r)353(p:^ - p„) 

X (qakaU-m„-^^j|qa), (13) 

where Ta{z) is the two-body scattering amphtude, and 

maimg + m^) mam^ ,^ ,, 

M = mo, + m^ + m^ , M„ = °V /J ^^ ^ ^^ = ^^_ ^ 14 

nia + mj^ + m^ mp + m^ 

we finally obtain 

(pLqa|m„/3(2:)|p/3q/3) = 5„/3(27r)^5^(p'^ - p„)(q^|T„( 2 - m„ - ^]^ J IQ^) 



Z.^^ ^a,jy ^27r)3 (27r)3 (27r)3 (27r)3 ^^""^ ^ ^^^ PJ^QJ^"!^^ '^^ 2M, 



(2vr)3^3(p. _ p^^(p:;,q:;))(2vr)3^3(q. _ q;;(p:;, g^;)) ^^ 

X 7-77T2 TZiFvi {V^^.^\T^ip[z)\ppi^p) . (15) 

M + ^^^ + ^^i^ - ^ - iO 
2M„ 2/i„ 

The relations between the momenta Pa,P/3,P7 and qa,q/3,q7 for different channels are given 
by 



nia rrijiirna + mp + m^) m^ 

\ P7 "•" ^1 1 *la — ~7 ; TT ; r P7 ; 



rria m^{ma + mp + m^) mp 

Pa = P/3-q/3, Qa = ^^ 77 ^ ^ P/3 ^ q/3 • (16) 

ma + m^ {ma + m^)[mi3 + m^) mj^ + m^ 

Here, (a/37) = (123), (231), (312). 

Further, one may express the effective one-channel potential through the solution of Faddeev 
equations. Removing first the CM motion, one gets 

(k;k'2|W(^)|kik2) = {2Txf5\k!, + k'2 - ki - k2)(qX;^)|q3) , (17) 



with 



3 

(q'|w(^)|q) = -y3(q',q)+ Y. J 



(27r)3 (27r)3 (27r)3 (27r)3 
/ II ii\i II ii\„ I \\ III iii\ ff^ I III iii\\* 



^J (27r)3 (27r)3 (27r)3 (2vr)3 

T' I I II ll\l II ll\„ I \\ III lll\ /T^ ( III lll\\* 

In the above expressions, the following notations are used 



ra(q;Paqa) = r q, -q; kc, = Pa, k^ = qa ■ p^,, k^ = -q^ ^ Pa i , 

\ m^ + m^ mp + m^ ' 

^a(q',q) = K.(q',-q';q,-q), (19) 

and 

iiiin^ (,MW'W"\-^ A (2vrf^^(pg-p;(0(2vr)3^3(q._q..) 

\Paqalg3,a/3l^j|P/3q/3/ " (>a3(>f33 7-77^2 rZ7i\2 

2M3 2/i3 

+ /„//N2 77/^2 (Paqalma/^l^jIP/^q/?) / ///\2 /^///n2 (^O) 

2Ma 2/ia 2M;3 2/i;3 

Finally, the LS equation for the two-body scattering matrix Tp{z) after removing the CM 
motion 

(k;k'2|Tp(^)|kik2) = (27r)353(k; + k'2 - ki - k2)(q:,|tp(^)|q3) , (21) 

takes the form 

/ /, / M ^ , /, / M V MV (q'|w(^)|q")(q"|tp(z)|q) , , 

(q'|tp(^)|q) = (q'|w(z)|q) + ^ 7^ q^ )^ ' ^^ 

mi + 1712 H z — iO 

2/i3 

The above formulas simplify considerably, if the pair potentials have separable form. For 
completeness, in Appendix A we list the pertinent expressions in the separable model for the 
case of three identical particles. 

2.2 The model in a finite volume 

Now, let us put the system described by the Hamiltonian in eq. (1), in a finite cubic box of a 
size L. Assuming periodic boundary conditions, the momenta of all free particles take discrete 
values ka = 27rna/L , a = 1, 2, 3 and ria G 1^^. The only difference between the infinite-volume 
and finite-volume cases consists in replacing the momentum-space integrals in all scattering 
equations by the sums over the discrete momenta. 

7 



The finite-volume counterpart of tlie LS equation witli tlie effective two-body potential in 
eq. (22) is given by 

(q'|tfWlq)Hq'|w^(.)|q> + i5: Wl^^f-'IO '^f^^^^'^ . (23) 

q" mi +1712-1 z 

Here and below, we attach the superscript "L" to the quantities defined in a finite volume. 
Note also that in the CM frame, the relative three momentum in the intermediate state is given 
by k'/ = — k2 = q". This means, that the summation momentum q" takes the discrete values 
q" = 27rn/L, n G Z^ 

The scattering amplitude tp{z) has an infinite tower of poles, corresponding to the discrete 
energy spectrum of a system in a finite box. It can be shown that the locations of these poles are 
determined by the Liischer formula. Indeed, performing the partial-wave expansion in eq. (23) 



(q'|t^(2;)|q) = 47r ^ Yi,^,{ci')tji^,^i^{q' ,q; z)Yil{q) , 

I'm' ,lm 

(q'|w^(2)|q) = 47r ^ Yi,„,,{q')w^^,^i^{q' ,q- z)Y;^{q) , (24) 

I'm' ,lm 

where q denotes the unit vector in the direction of q and l^m(q) stands for the spherical 
function. The position of the poles is determined by the equation 

detl) = 0, 

l^l'm'M = 5vi5m'm-y^^^:—Kl;^^,iu^u{p,p]z{p))Ml'im".lm{l^), 

^ — ' ZTT ' 

l"m" 

(-)'' ''"^' ^ V 

■Ml'm',lmi^) = 3/2 / , / ^ ~~n"^i'^(-*-' ^ )Clim' ,js,lm ■ (25) 



^3/2 ^^ ^^ ^j- 
j=\l'~l\ s=^j 



Here, 



^ = — , (26) 

27r ' ^ ' 



the quantity Zjs{l] z/^) denotes the Liischer zeta-function, and the symbols Ci>m',js,im are given 
by (see, e.g. [27]) 

C,m',sM = (-)™/-^+V(2/' + l)(2j + 1)(2/ + 1) f£, i Mfj! i ^V (27) 



m' s —m / \ / ' 



Further, the quantity K^(q, q, z(q)) with z(q) = rrii + m2 H is the on-shell solution of the 

2/U3 



LS equation up to the exponentially suppressed contributions in the box size L. It is given by 



Kvm'MiQ^ ?; z{q)) = w^^,i^{q', q; z{q)) + -^ P.V. / 

27r Jq 



2 



Wr-q 

X Yl ^^m',l"n,"{(l'. Q"; 4Q))Kv'm",lm{(l". Q', z{q)) • (28) 

l"m" 

Here, P.V. stands for the principal- value integral. 

In the above equations, the mixing of the partial waves occurs, because the rotational 
symmetry is broken on the cubic lattice. The equation can be partially diagonalized in the 
basis of the irreducible representations of the cubic group [27,37]. In case of the fermions, 
the basis of the irreducible representations of the double cover of the cubic group should be 
used [38]. 

The expressions simplify considerably, if, e.g., one assumes that only the S-wave scattering 
contributes. Then, the Liischer formula can be rewritten as 

.3/2, 



jy) = -S^{p) + 7rn, n = 0,l,---, tan0(z/) = -^^^— ^ . (29) 



^oo(l; i^ 
The quantity S^{p) stands for the so-called pseudophase [30,32,33] 



tan5^(p) = ^<,oo(P,P;^(p))- (30) 

If the energy z is below the three-particle threshold, the volume-dependence in the effective 
two-body potential is exponentially suppressed in L. Consequently, up to the exponentially 
suppressed contributions, the pseudophase S^{p) does not depend on L and on the level index 
n (the latter dependence arises through the dependence of the effective potential on L = Ln{p) 
at a given value of p). In this case, 6^{p) = 6{p) coincides with the conventional elastic 
scattering phase. 

Above the 3-body threshold, the effective potential is given by (cf. eq. (18)) 



(q'|w^(z)|q) = -V3(q',q) + T^ J2 E ^ .(q'; pX 



L 

X (pXlg3,a,(^)|pW)(f/3(q;p"q"))*- (31) 

Note that in the above equation, we do not attach the superscript L to the quantities V3 and F, 
which are the same in a finite and in the infinite volumes. Further, the summation in eq. (31) 
runs over the momenta (cf. eq. (12)) 

27rn'' „, 27rn' 



,, ^""q 


p;^'=- 


■'"P 








p„ ^ ■ 


L ' 




q" = '" (y + 


771/3 + "1-y 




q^ 


271 

L 


fyu , "^7 




n" W" 1" V" ^ 


Z3. 











-M' 



(32) 



The finite-volume version of eq. (20) reads 

^ (P3L + iS3L_^ 
2M3 2/i3 

1 1 



(P::q'a|m^^(z)|p;q;') . „n2 ,^.n2 ' (33) 



2M„ 2/i„ ^ 2M^ ^ 2/i^ 

Faddeev equations in a finite volume are written as (cf. eq. (15)) 

1 ^ 

(Paqa|m^/3(^)|P/3q/3) = -^a/^^^^^p^ (q^ |t^(^; Pa) iQa) + "^ Xl^^ ~ ^"^) XI ^^V.: 

7=1 p^,q'^,p^',q;' 



X (qaka(2; Pa)|qa) ^^7W . /M2 (P.yq7lK/3(^)IP/3q/3) ' (34) 

^^ iPoL iqa^_^ 

2M„ 2/i„ 
where the two-body scattering amplitude in a finite volume is a solution of the equation 

/ / 1 L/ ^l \ ^ t7^^' ^ ^^^ ^ ^:^ (-'^a(qL^q"))(qal'r«(^;Pa)|q») /o^-^ 

(qaFa(2;;Pa)|qa) = (-^(qa, q^)) + "p 2^ -"2 f-jr^ ' (35) 

^ q^ M + -^ + ^^^ - 2 

2M« 2fla 

Note that the momentum p^ enters in the definition of r^ twice: in the argument z — )■ 
2; — rriry ^, and through the shifting of the summation over the momenta a'' which, 

2Mo, to to Ma > 

according to eq. (32), is no longer proportional to an integer number. In the infinite volume, the 
summation is replaced by an integration over the whole momentum space, and the dependence 
on Pa remains only in the argument. 

As in the infinite volume, the equations simplify considerably, if we assume pair interactions 
of the separable form. The pertinent equations are listed in Appendix B. 

2.3 Singularities of the effective two-body potential 

As mentioned above, below the three-particle threshold the finite- volume contributions to w^{z) 
are exponentially suppressed. Therefore, up to such suppressed terms, the effective potential 
coincides with the regular function "w{z), which is defined in the infinite volume. 

Above the three-particle threshold, the singularities emerge in w^{z). In the vicinity of 
these singularities the pseudophase rapidly changes by tt. These singularities may strongly 
affect the finite volume spectrum in the vicinity and above the inelastic threshold. For this 
reason, below we shall discuss them in detail. 

Potentially, the matrix element (q'|w^(2;)|q) may become singular, when 

10 



E ^St + Z 



+ /. ^v — + z^ a:v/c?v + 



a opz^p 




Figure 1: A multiple-scattering series for the three-body Green's function. 

i) the three-particle denominators M -\ -^ — I 2; in eq. (33) vanish; 

ii) the matrix elements of the operator m^^(2;), which are solutions of the Faddeev equations 
in a finite volume, develop a pole at those values of energy, which do not coincide with 
the poles of the three-particle energy denominators. 

We are going to demonstrate below that the singularities of the first type cancel, and only of 
the second type survive. In order to prove this statement, we mention that, from eqs. (32) and 
(35) it follows that the two-body scattering matrix t^{z) vanishes exactly for those values of z 
where the three-particle propagators develop a pole. These energies are given by the following 
equation 

^ = M+^ + -^fp; + ^^p'^y, (p:,p:;) = ^«,n:;), n::,n:;6Z^ (36) 
2Ma 2fia \ mp + 777,7 / ' L ' ' 

Let us show now that there are no three-particle singularities of the first type in the matrix 
elements of the three-particle Green's function. For simplicity, we check this property only for 
the lowest three-body singularity aX z = M. 

Let us consider the multiple scattering series for the three-body Green's function which is 
shown in Fig. 1. According to this, the sum of the matrix elements of gj^^siz) is written as (cf. 
eq. (33)) 

3 r3x , 7-3/: , 



a,/3=l M + -^ + ^^ 



2M3 2fl. 



3 






2M3 2/i3 J\ 2M3 2/i3 

Here, we have used the property 



11 



The first term in eq. (37) has a pole at z = M, when P3 = Ps = qs = qs = 0. The denominator 
in the second term has a double zero here, but the matrix element of the operator t^{z; p^) 
also has a zero, so that only a single pole remains. This property can be easily checked in 
higher orders. Introducing the notation 



^""L3 M-z ' ^'^^^ 



we obtain (see eqs. (33) and (34)) 



^ (00|g3^„^(2;)|00) = -I^^!^l + (y,+y, + y,) 



a,l3=l 

+ (2/1(^2 + 2/3) + y2{yi + y-i) + ysiyi + 2/2)) + ■ ■ ■ } + non-sing , (40) 

where "non-sing" stands for the terms that are non-singular at z = M. 
Summing up the series (see Appendix C), we find 



3 



L^ XiX-zXj, 



J2 mgL,i^)m = - "7^"^ . + non-sing, (41) 

^, M-z a;iX2 + X2a;3 + 0:3X1 - 2x1X2X3 

a,fi=l 

where Xa = 1 + ya- Let us now show that, in fact, the first term in eq. (41) is also non-singular. 
To this end, we consider the equation for the operator r^ 

(q>i(^. 0)M = (-K«, qj) + ^ E '-^°'<-<>;<y-^°>l'^°> . (42) 

^ q^ M + ^^^ - Z 

In order to separate the singularity at ^ = M, we single out the term q^ = in the above 
sum. It is then straightforward to see that 

where 

(q>L(^;0)lq„) = (-V-.^.q,)) + ^^ ^ <'-^'^'i^'<))\'f>'^^<'-''')\'^-) ^ (44) 



2/i„ 



Since the term with q^ = is absent in eq. (44), the quantity {c(^\T^^{z]Qi)\(\a) is non- 
singular at 2 — ;■ M. Now, the validity of our statement after eq. (41) follows from the fact that 
Xa = '\- + ya = 0{z — M) as z — )■ M and, consequently, the three-particle singularities cancel in 
the sum of the matrix elements of g,^^p{z). 

In conclusion of this section we summarize the main points: 

12 



i) The singularities of the effective two-body potential y^^{z) are caused by the singularities 
of the matrix elements of the operator \\\aji{z) which emerge as a result of solution of 
Faddeev equations in a finite volume and thus have a non-perturbative origin. 

ii) In other words, the singularities do not emerge, when the energy coincides with the eigen- 
values of the free Hamiltonian in the box. The singularity structure is determined by the 
eigenvalues of the full Hamiltonian. 

iii) Individual terms in the finite-volume multiple-scattering series, however, contain the sin- 
gularities determined by the free Hamiltonian. In order to make them disappear, multiple- 
scattering series should be summed up to all orders. The Liischer formula ensures such a 
summation in the two-particle case. In this section we have explicitly demonstrated that 
the summation can cure the problem in the three-particle case as well. 

iv) The results of this section may serve as a warning against an approximate truncation of 
the multiple-scattering series in a finite volume. While in the infinite volume this may 
affect only the numerical precision, the singularity structure of the effective potential can 
be modified drastically as a result of such a truncation in a finite volume. Consequently, 
if the approximations are made in the three-body case, the singularity structure of the 
effective potential should be examined very carefully. 

3 An alternative derivation of the Liischer formula 

In order to find the energy spectrum in a finite volume, the equations, which are considered 
in section 2.2, can be solved numerically, applying, e.g., the method of ref. [39]. However, 
these equations contain potentials as well as off-shell two-body scattering matrices, which are 
model-dependent. The central question is, whether the predicted energy levels are also model- 
dependent. In other words, if two different potential models lead to the same S'-matrix in the 
infinite volume, can the finite-volume spectra in these models be different? 

In case of the two-particle elastic scattering, the answer is given by the Liischer formula, 
which relates the finite- volume spectrum to the (on-shell) S"- matrix element. Our aim is to 
rewrite the three-particle equations in a finite volume in a similar fashion, in terms of the on- 
shell ^'-matrix elements only. In order to do this, in this section we consider a novel derivation 
of the Liischer formula. The method used here can be generalized for the case of three particles, 
as shown in section 4. 

Let us consider the sum 

where $(p) denotes a regular function of p. Next, we perform a partial-wave expansion 



Im 






13 



where yimip) = p'^«m(p)- The sum in eq. (45) can be rewritten in the following form 

^0 ^ p im P ^0 

Note that we have introduced the momentum cutoff A, in order to regularize intermediate 
expressions. The cutoff disappears from the final expressions^, if the function $(p) falls off suf- 
ficiently fast with p Further, we have introduced a regulator ][x) with the following properties: 

1. The function /(x) is bounded and smooth (has any number of derivatives) on the whole 
interval x G ] — oo, oo[. 

2. f{x) = 1 if x> 0. 

3. lima._j.o- /(3;) = 1 and lim2._;.o- /'■"^(x) = for all n 7^ 0. 

4. The function j[x) vanishes exponentially when x — )■ — 00. 
Otherwise, the function /(x) is arbitrary. An example of such a function is 



f{^)={ ..J 1 ^ .. ./n (4^ 



1 


1 \ 




x> 
x<0 


cxp 


- exp(x"2)y 



The choice of the scale /i in the regulator / is also arbitrary. For example, one could choose [i 
to coincide with the mass of the lightest particle. 

Taking now into account the fact that the expression $;(p)/p' is a regular function of p^ 
in the vicinity of p^ = 0, the regular summation theorem [40] can be applied. Up to the 
exponentially suppressed terms in L, one may replace the first sum in eq. (47) by the integral 

IpKA rf3p 3;^^(p) /$^(p) $^(^^) ^^ 



' ~ J (^trp'-^o-^ov^^ ~ir^^''^^ 

+ ^/(.oV/^^)7^'EE¥4, (49) 

^0 ^ p tr p ^0 

where Qq -^ Qo + iO prescription has been chosen arbitrarily (the numerator of the integrand 
vanishes at p"^ = q^, so the prescription does not matter). Simplifying the above expression, we 



^At finite values of A, rapidly oscillating terms at i — > 00 may occur for the sharp cutoff, so a mathematically 
rigorous procedure is to use a smooth cutoff at a momentum scale A. In order to ease the notations, we however 
proceed further with s sharp cutoff. The oscillating terms are briefly considered in Appendix D. It is shown 
there that these - as expected - are harmless. 
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arrive at 






^/(%V^^)E^^'-(1'-')' (50) 



47r2L 

Im 

where u = qQL/{27T), and Zim{l] u"^) stands for the Liischer zeta-function 



Zir^il; u') = hm ( V ^(A^ - n^) ^^^ - dio^mo^AnX \, A = ^ . (51) 



■^^ = hm ^ V ^(A^ - n^) ^^ - 5.0^^0 v^AJ , A = ^ 

A^oD -^^^ n^ — z/^ 27r 



Formally, the above equation can be rewritten in the following manner 

^^(2^pj^ = p.v.-^ + 5:-i,y,;(p)^,„.(l;.Xp^,^. 

P'V-^ - ;;i4^-"^(P^«S. (52) 

where the quantity A(p^, ^q) coincides with the conventional Dirac (5-function 5(p^ — ^q) for 
Qq > and is defined through the action on the smooth test functions 

--/ (^ A(p^go^)$(p) = ^^^^%{%)Ml^^') (53) 

(recall that /(^'o//^^) = 1 for gg > 0). For q^ < 0, the action of the distribution A on a test 
function is defined through the analytic continuation of $o(q'o) in eq. (53). The regulator / 
serves the purpose to effectively cut the contributions with —q^ > fi"^. 

To summarize, the momentum-space two-body Green's function in a finite volume, which 
is given by 

«j,k;.).|5:M;|!(-jE), ,54) 

-'^ p p % 

where u is the reduced mass and z = mi + m-i H , can be decomposed into the following 

2/x 

^7rA(k^ ql) + A(k^ ql) J^ F,;.(k)-^ Z,„(l; u') 

Im 

= Go(k;^) + Gu(k;^) + GF(k;2) = GK(k;2) + GF(k;^). (55) 

The equation (55) gives the splitting of the finite-volume two-body Green's function into the 
infinite-volume part Gk = Gq + G\j (corresponding to the principal-value prescription at the 
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;s 








Gj(kiz) = 


= *|k^- 


1 


-zO 



singularity) and the correction Gp- The latter is proportional to the on-shell factor A(k^,gQ), 
which projects onto the energy shell. This fact plays a crucial role in the proof of the statement 
that the finite-volume spectrum is determined only by the infinite-volume S'-matrix elements 
which are defined on shell. 

The Liischer formula can be derived straightforwardly, using eq. (55). To this end, we 
take gg > 0. Then, A(k^,gQ) is replaced by the conventional (5-function. First, note that, 
substituting the free Green's function by the quantity defined in eq. (54), the two-body LS 
equation in a finite volume can be written formally in the same form as in the infinite volume 

(p|T^(z)|q) = (p|(-V)|q) + y ^^(p|(-V)|k)G^(k;z)(k|T^(z)|q). (56) 

Next, we define 

T = (-V) + (-V)GoT, 

K = (-V) + (-V)(Go + Gu)K = T + TGuK, 

T^ = (-V) + (-V)(Go + Gu + Gf)T^ = K + KGfT^. (57) 

It is immediately seen that T and K are the two-body T- and A'-matrices in the infinite volume, 
respectively. The relation between the /f-matrix in the infinite volume and the T-matrix in a 
finite volume T^ is given by the last equation in eq. (57). After carrying out the partial- wave 
expansion 

(p|T^(z)|q) = 47r5^5^y,^,(p)T,f„,,,^(p,g;z)l^;(q), 

I'm' Im 

(p|K(^)|q) = A7rJ2yim{p)Ki{p,q;z)Y;^{ci), (58) 

Im 

and integrating over the angles, this equation can be rewritten in the algebraic form 

Tvm'M^P^ 1'^ ^) = 6w6mm'Kl{p, 9' ^) + ^ ^'''(P' ^O; z) ^ A^i'm',i"m" ('^)T;^m",im(90, 9, ^o) ' 

l"m" 

(59) 

where the matrix Aii'm',imi^) ^^ displayed in eq. (25). Going to the mass shell p = q = qo 
and requiring that the above system of linear equations is singular {T^ becomes infinity, that 
means that the determinant of the above system of linear equations vanishes), we finally arrive 
at the Liischer formula, see eq. (25). Note that, since in case of an elastic two-body scattering 
considered here, the i^-matrix in the Liischer formula stands for the infinite- volume i^-matrix, 
^Vm' Im ^^ ^1- (25) should be substituted by SwSmm'Ki. Note also that, within the normalization 
used, the on-shell /T- matrix is related to the scattering phase shift, according to 

LiT) 

tSin6i{qo) = --Ki{qo,qo]z). (60) 

zir 

In conclusion, several remarks are in order: 
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The decomposition of the two-particle Green's function in a finite volume, given in eq. (55), 
is valid, if both sides of the same equation are integrated with the same regular test 
function. In general, it does not hold, if integrated with singular functions. In the Born 
series of the LS equation, the two-particle Green's function is always integrated with the 
potential, whose singularities are assumed to lie far away from the integration contour. 



ii) If v^ < 0, then Zimil; v^) = —SioSmo'^^^'^V—i^^ — iO, up to the exponentially suppressed 
contributions. This implies that the sum of all the terms, proportional to $o(q'o) for i^^ < 
in eq. (50), is in fact exponentially suppressed. Hence, the quantity 5*2 does not depend 
on the values of the function $(p) outside the original range |p| > 0, as it should. 

iii) The reason why we still retain these terms and do not replace the quantity A(p^,gQ) by 
(5(p^ — Qq) everywhere, is the fact that the principal-value integral, which was defined 
above, is the regular function of Qq (no unitary cusp present). This property plays no role 
in the two-particle scattering (even in the multi-channel case) but, as we shall see below, 
becomes critical in case of three-particles. Note also that in ref . [32] , which deals with the 
multi-channel scattering for the coupled -ktt-KK channels, we have effectively used the 
same splitting as in eq. (55) but with f{x) = 1. 

4 Three particles in a finite volume 

4.1 Splitting of the three-particle Green's function 

In analogy with eq. (54) we define the finite-volume three-body Green's function in a channel 
q; = 1, 2, 3 as (in order to ease the notations, we suppress the channel index a in all momenta) 



Goa(k,i;^) = ^E 



(27r)353(p-k)(27r)353(q-l) 



L"" ^ ,, p' q' 

rriR 2n _ _ ^q 

q = P + P, (P,P) = ^(n,n), n,nGZ% (61) 

m/3 -|- ?n^ L 

so that, in the above equation, the summation is carried out over the integers n, n. In order to 
find the desired form of the splitting, we separate the 2-particle Green's function in the above 
expression 

p q ^ua 

< = 2/^-(^-^-^)- (62) 
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Acting in the similar way as in the two-body case (see section 3), we obtain 

1 



GL{K^;z) = ;^ E(2^)'^'(P - k) { 



M + -^ + z-iO 



2M„ 2/ic 



«7rA(l2,2/i„(M+ ^ 



2M„ 



+ A(^1^2;x„(M+^-^))5^r,:;,(i)Az^Jl;z.^)}, (63) 

where z/^ = qoaL/{27i), a = [niis/^mis + m^)] pL/(27r), and the Liischer zeta-function in the 
moving frame is defined as (cf. the pertinent expression in the CM frame, eq. (51)) 

^L(l; ^^) = ^li^m I J2 0i^" - (n + a)^) .^'"/^.^-l^ - ^^o^^oV^aJ . (64) 

Since the principal-value integral over the variable 1, containing a regular function of the 
arguments p, 1, is a regular function of the remaining variable p, one may use the regular 
summation theorem in this variable and rewrite eq. (63) in the following form (cf. eq. (55)) 

GUk,\;z) = p ^2 ^^A(l2,2/iJM + A_-2;)) 

M+-^ + — -z-zO ^ ^ ^^^^" ^^ 

2M„ 2fia 



= Goa(k, 1; z) + GuaiK 1; z) + GFa(k, 1; z) 

= GKa{Kl;z) + ^J2^27rf^'{p-k)G,4p,\;z). (65) 

p 

Needless to say that the above splitting is valid, if both sides of eq. (65) are integrated with 
a regular function of the momentum variables k and 1. Further, since the variables k and p 
are not restricted from above, the argument of the distribution A can become positive and 
arbitrarily large, independent of the choice of the variables P and z. This means that one has 
to necessarily deal with the analytic continuation of the pertinent amplitudes below threshold. 
Note however that this sub-threshold contribution is effectively cut by the regulator / for the 
momenta —qla > /^^- Assuming the range of the potentials much smaller than the inverse of the 
lightest mass in the system, one does not expect to encounter any singularities in the analytic 
continuation for —q^^ < /i^. The above argumentation serves to justify the introduction of the 
regulator /. 

Using the representation of the three-particle Green's function, given in eq. (65), one may 
try to repeat the same steps as in eq. (57), also for the three-particle T-matrix. In the three- 
particle case, however, a new complication arises, related to the presence of the disconnected 
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contributions (the diagrams where a spectator particle propagates freely, when the other two 
particles interact). These diagrams contain the factor L^Sp^^^ ({27rY6^{pa — Qa)) in a finite 
(infinite) volume. This factor is not a regular function and is L-dependent (in a finite volume). 
Consequently, the method, which was used in the two-particle case, can be applied here first 
after these disconnected contributions are removed. Below we shall demonstrate, how this can 
be achieved. 

4.2 Effective two-body A'-niatrix in the presence of the 
three-particle intermediate states 

As was shown in section 3, in order to obtain the Liischer formula from eq. (23), one may 
substitute here the splitting of the two-particle Green's function, given in eq. (55). The result 
is given in eq. (25). The effective two-body if- matrix obeys the equation 



(,lK^(.)|.)H.1w',.,|<,)H-P.v7ig<3VM™S!(£M. (,,) 

mi +1712 + -^ z 

2/i3 

Performing the partial-wave expansion of the above equation, we arrive at eq. (28). Note that, 
in difference to eq. (59), which refers to the elastic case, the effective i^-matrix in the above 
equation still depends on L above the three-particle threshold. 

Taking into account the expression for the potential w^{z), given in eq. (31), we get 

(q'|K^(.)|q) = (q'|K3(^)|q) + -^ Yl r„(q';pX) 

n"n"n"'n"' 

X (P::q::ig3"(^)|p"q")(f/3(q;p;^'q;^'))% (67) 



where K^ is defined in the infinite volume 



rfV {-V^{q',q")){q"\K,{z)\q) 



(q'|K3(.)|q) = (-l^3(q',q)) + P.V. / 7^ ^^^^^^^^^^^^^^^TTj^^^^ • (68) 

7711 + 7712 + — Z 

2/i3 

Further, the quantity F^ is defined through F^ in the following manner 

r ( ' "r."\ V ( ' "r."\ ^ f ^'^" (qlK3(z)|q^^)r.(q^pX) ..Q^ 

ra(q ; P^qa) = r^(q ; P^q^) + / 7^^ 7-77T2 > (69) 

mi + m2 + — z 

2/i3 

and the Green's function g^iz) obeys the following equation 

(Paqalg3(2^)|p/3q/3) = (Paqalgf(2;)|p/3q/3) + 7i^ Yl (Paqalg^(2)|P7q7) 

x(p:;q:;i(-v4(.))ipX')(pXlg3(-~)|p/5q/3) • (7o) 
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Here, g^i^z) = X^a /3=i Ssa/^l'^) ^^^ 7'"^ "^^^ t)e chosen arbitrarily. Further, 'V4,{z) denotes the 
effective three-particle potential, which has emerged after "projecting out" the two-particle 
intermediate state. This potential is given by the following expression 



(q'|K3(z)|q) 



mi -F m2 -h — 2 mi -F 7712 -h 



2/i3 / V 2/i3 



q2 

7711+7712 + Z 

2/i3 



r/3(q;P/3q/3)- (71) 



We would like to stress that the effective potential y^A{z) is defined in the infinite volume and 
is a smooth function of its arguments. 

To summarize, after "projecting out" the two-particle intermediate state, an effective three- 
particle force V4(z) appears in Faddeev equations, even though we initially assumed that no 
such force is present. The effective two-body i^-matrix in a finite volume is given explicitly 
by eq. (67). In order to evaluate the finite- volume effects in the three-body Green's function, 
which enters this expression, one should consider Faddeev equations in a finite volume. 

4.3 The three-particle counterpart of the Liischer formula 

In order to derive the three-particle analog of the Liischer formula, we have to deal with the 
three-body equations in the presence of the three-particle force. To this end, one may use, 
e.g., the formalism described in the papers [41] - with the potentials replaced by the ii'- matrix 
elements and the free Green's function replaced by Gp. However, as discussed above, there 
exists a problem related to the presence of the disconnected parts. Namely, in the presence of 
the Kronecker-(5, contained in the disconnected parts, the use of the splitting procedure for the 
three-particle propagator according to eq. (65) can not be justified mathematically. In order to 
circumvent this problem, we act by using a trial and error method. Namely, we first apply the 
splitting in the three-body LS equations, as if there were no disconnected parts, and further 
use the Faddeev trick. In the resulting equations, the disconnected terms, containing the 6- 
functions, emerge in a finite volume. At the next stage, we discard these singular terms by hand, 
thus making a conjecture about the correct form of the equations. At the final step, we check 
this conjecture explicitly, by considering the multiple-scattering series that emerge from the 
resulting equations, and showing that this series coincides with the original multiple-scattering 
series in a finite volume. 

Symbolically, the result for the effective two-body i^-matrix can be written as follows 

K^ = K2-,2 + K2^3(Gf + GfRfGf)K3^2 , (72) 

where Kj_^j denote the pertinent i^- matrix elements in the infinite volume'^, Gf stands for the 
finite-volume part of the three-particle Green's function (see eq. (65)), and the quantity Rf is 



^ Above threshold, the three-body X-matrix coincides with the one defined, e.g., in refs. [42,43]. 
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given by 



4 3 

Rf = 2_^ ^f^i^ + z^ ^° 

fi,iy=l a=l 



R4/3 — O^Gf ( ^/3 + 2_^ -^7/3 j ' 
\ 7=1 ^ 

R-Q/S = ^aGp ( ^^(1 ^ ^a7)R7/3 + R4/3 1 5 
^7=1 ^ 

3 
Rq,4 ^ C/q,(jif / (,1 — 0a7jR74 ~l~ C'aVJFR44 



7=1 



3 

R44 = ^4 + 04Gf y R^4 . 

7=1 






Here, /i, i/ = 1, ■ ■ ■ ,4, whereas a, /3, 7 = 1, ■ ■ ■ , 3 and 






0^ = K^ + K^Gf^^ , K^ = (-V^) + (-V^)(Go + Gu)K^ , 


K3_>3 -- 


4 



(73) 



(74) 

Note that in eq. (73) (third hne), we have omitted the terms of the type ^qGf^/? with a ^ p. 
Physically, such terms correspond to the finite-volume corrections in the disconnected diagrams 
and emerge, if one faithfully applies the splitting procedure even to the disconnected piece. This 
omission will be justified below by showing that eq. (73) produces - diagram by diagram - the 
correct splitting of the infinite- and finite- volume parts in the multiple-scattering series'^. Here 
we also would like to mention that this system of equations is formally identical to the equations 
that relate the T-matrix and /f-matrix elements, with a replacement R — )■ K, K — )■ T and 
Gf^Gu. 

The above equations can be written down in the explicit form 

(q'|K^|q) = (q |K2-.2|q) + (q'|K2^3GFK3->2|q) + (q'|K2^3GFRFGFK3^2|q) , (75) 

where^ 

(q'|K2-^3GFK3^2|q) = J^^ 7^^ {(^\'i^i%\Pa(la)GFa{Pa, qa] z){pa(la\K.^^J,2h) > (76) 

Pa 



^The physical meaning of this prescription is very transparent. The finite- volume corrections emerge only 
in the loop diagrams. However, due to the presence of the Kronecker-delta in the disconnected diagrams, a 
first iteration of the disconnected diagrams in the Faddeev equations gives a connected diagram without a 
loop. Its explicit expression is identical in a finite and the infinite volumes. The loops (and, consequently, the 
finite- volume corrections) emerge first in the second iteration, see fig. 2 

""Note the superscript (a) in the i^-matrix elements. It emerges because Guq explicitly depends on the 
channel a through the regulator f{qQa/^J■^), see eq. (65). 
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Figure 2: a) One iteration of the disconnected diagrams produces a connected diagram with 
no loops; b) Loops emerge first after two iterations. 

The last term in eq. (75) is given by 

(q'|K2^3GFRFGFK3^2|q) = (qlK2^3GF('X^0«+ X^R^.'jGFKg^alq), (77) 



where 



(q'|K2^3GF0aGFK3^2|q) = ^J2j 



(27r)3 (27r)3 ^^ l^243lPaqa)G'Fa(Pa, q^! A 



(a) 



and 



X (qal^"(^;P")|qa)G'Fa(pa,qa;2)(paqa|Kn2lq) 



(78) 



(q'|K2^3GFR,.GFK3^2|q) = ;^ E / \^Y in)3 (qlK^l3lpU;)gF.(pL, q:.; z 



L^ ^ J (27r)3 (27r)3 



(/3) 



X (Paqa|RA'^|P/3q/3)G'F/3(P/3, ^.p] z)(p/3q/3|K^42lq) 



(79) 



Note that if /i = 4 and/or i^ = 4 above, the pertinent indexes a, /3 can be chosen arbitrarily. 
Next, 6a are determined through 



(qal^a(^;p«)|qQ) = (qU^af ^^-"^a-^^jlqa 



rf3l^ 



q'„|Kc, z-nia 



_pL 

2M„ 



|la)G'Fa(Pa,la;2)(la|^a(2;;Pa)|qa) 



(80) 



In this expression, Kq stands for the two-particle Ji'-matrix. Performing partial- wave expansion 
in this equation and integrating over the absolute value |1q,|, we arrive at a set of algebraic 
equations that relate Kq, and 6 a- These equations mix all partial waves. Hence, in order to 
solve them, a partial- wave truncation is necessary. It is easy to recognize that eq. (80) is nothing 
but the Liischer formula for a two-particle sub-system in the moving frame (the third particle 
plays the role of a spectator). The result was of course expected from the beginning. 
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The equation for ^4 is given by 

(Paqal^4|P/3q/3) = (p^q'^|K4|p/3q/3) 



1 . r d'^i 

+ -j^ / TT-^ (p'aq'„|K4|k„la)GFQ(k„,l„;2)(k„l„|04|p/3q/3). (81) 



It is easy to see that, after performing the integration over the absolute value of the relative 
momentum !„, eq. (81) transforms into the matrix equation that can be solved straightfor- 
wardly by using numerical methods. The matrix indexes here include the discrete values of the 
spectator momenta p^, q/3, ko,. In order to solve this equation, a method analogous to the one 
described in ref. [39] can be used. 

Finally, we write down the equations (73) in the explicit form 

(Paqa|R'4/3|P/3q/3) = / —^ {p'M^i\Pl3(l'li)GF(3{Pl3, q!^] z){q!^\0(s{z;P(s)\q,3) 



(p;q:,|R„;3|P/3q/3) = X^(l - 6^,) I ^ K\e{p'a, z)\c^:,)GUp'a. <; ^)(p:;q:;|R,/3|P/3q/3) 

+ l^,Meip'^-z)mGUp'a,<-,^){p'a<\^^^\Ppqp), (83) 

(Paqa|Ra4|P/3q/3) = ^(1 - (5„^) / — ^ (q^ |6>(p^; z) |q;^)GFa(p^, q^; 2) (p:;q:j|R^4|p^q;3) 

+ y"^(q:,|e(pL;^)|q::)GF.(p:„q'^;^)(pX|R44|P/3q/3), (84) 



(pLqa|R44|P/3q/3) = (Paqal^44|p^q/3) 



1 f da 

+ EtsE / 7^(Paqal^44|pX)'^F7(Prqr^)(pXI^74|P/3q/3) 

7=1 p;^' ^ '"^^ 



i5) 



where 



P7 = : Pa-q^^ % = 7 \ \7 \ \Pa : qo- l^oj 
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4.4 The term without a potential insertion 

Below we shall start examining the multiple scattering series which are obtained from eq. (67) 
diagram by diagram, and shall demonstrate - for a few illuminating cases - that our equations 
indeed produce the desired splitting of the infinite- and finite-volume contributions. 

We start from the trivial case of a diagram without rescattering, shown in fig. la. This 
diagram should be folded with the vertex functions F, F*. As a result, we arrive at the following 
expression 

J 1 ^^^ r3(q^;p3q3)(f3(q;p3q3))* 

^0 = I^ Z^ ^2 -2 • (87) 

psqs M + -t^ + -^ - Z 

2M3 2/i3 

This expression is of the type already considered in section 4.1, since the product of two vertices 
is a regular function of the momenta. Consequently, up to exponentially suppressed terms, 

T = P\l f ^^P=^ ^^^3 r3(q^p3q3)(f3(q;p3q3))* 

^ 2M3 ^ 2/^3 
+ ;^E/^^3(q';P3q3)GF3(P3,q3;^)(f3(q;P3q3))*. (88) 



4.5 Disconnected contributions 

The corresponding diagram is shown in fig. lb. Expanding Tq in Born series, one obtains the 
diagrams with one, two, . . . insertions of the potential. The explicit expression of a diagram 
with one insertion is given by 

rd^X^ ]_ \^ r^(q'; Paq.'a){-Va{q'a, qa))(ra(q; Paq^))* , . 

V 2M^ 2/i„ J V 2M^ 2^a J 

Carrying out the splitting of the infinite- and finite-volume contributions explicitly, it is easy 
to see that, up to the exponentially suppressed terms, the above equation can be rewritten as 



^1 = Z^y (^^^^^r„(q;p„qJGKa(p„,q„;2:)(-K.(q„,qa)) 

X G'Ka(Pa,qa;2;)(fc,(q;Paqc,))* 

^1 f (Pq' cPq 

+ 5Z X3 XI / (-27r)"3 i2-KY ^"*^^'' P"^")^!"*^^'"' ^"' P°' ^)(^°(^5 Paqa))* , (90) 



where 

Uti^a^^c.;Pa,z) = GKa(Pa,qa;2;)(-K.(qa,qa))G'Fa(Pa,qa;2;) 
+ GfaiPa, q^; Z){-Va{q'a, (la))GKa{Pa, (la, z) 
+ G'Fa(Pa,q^;2)(-14(qa,qa))G'Fa(Pa,qa;2;) . (91) 

The generahzation for the diagrams with many potential insertions between the same two 
particles is straightforward. Further, recalling the definition of K2^3 and K3_5.2, it is easy to 
ensure that the diagrams considered in sections 4.4 and 4.5, can be unambiguously identified 
with the pertinent diagrams emerging in the perturbative expansion of the quantity K2-^2 + 

K2^3(Gf + Gf ELi ^aGF)K3^2. 

4.6 Connected contributions 

The connected contributions emerge as a result of the expansion of the diagrams shown in 
fig. Ic in powers of the potentials Va- The lowest-order term has O(V^) and is given by 

1 f „(q'; p'aO{-vMa, P/3 + V7^^^^ p:.)) 

2 Z^ 2^12 ^ / , "^/3 , x2 

V 2M, 2/i, ;V 2M„ 2/i„ 

— Tfi ~ 

(-^/3(-Pa - ,. P/3, q/3))(r/3(q; P/jq/?))* 
X .^''7 . . • (92) 

^ _p^ _q^_ 

2M;3 2/i^ 
Up to the exponentially suppressed terms, this expression can be rewritten as 
d^P'a d^^a d^Pp d^q^ 



^2 = E/|^^^^r„(q'P«q'«)f^-.(p:^q:^'P/^q/^)(r/^(q'P/^q/^))* 

+ Ei^E /^^^r.(q'PX)f^i(pX;P.q.)(f,(q;P,q,))* 

+ EiiE/l^^^ r„(q ' p:.q'J^J(p«q«; p.q/3)(f ,(q; p,q,))* 

+ E 76 E / ^^ r„(q'; pLqL)t/2f,(p'.q;; P^q/^)(f .(q; P,q,))* , (93) 

where 

^2i(PaqL; P/3q/3) = <^Ka(p'„, qL; 2)(-'^^a(qa. qa))GKaip'a, ^a, z){-Vp{q'^' , q/3))GK/3(P/3, q/3; z) , 

(94) 
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+ GpalPa, q^; 2)(-K«(qa, qa))G'Ka(pL, q'^; 2^) (-^^/^(q^', q/3))GK/3(p/3, q^] z) 

+ G'Fa(Pa, qa; z){-V^{ql^, c{^))G^o.{p'^, q^; z)(-\/^(q;^', q;3))GK/3(p/3, q/?; 2) , 

(95) 



^2a/3(Poq'a; P/^q/?) = G'Ka(Pa, q^; 2;)(-K(qa, qa))G'K/3(P/3, qj^'; 2:) (-1^/3 (q^', q/3))G'F^(p/3, q/3; z) 

+ GkcIpL, qL; 2;)(-V'a(qL, qa))<5F/3(p/3, q^'; 2)(-V"/3(q^', q;3))G'F/3(p/3, q/?; z), 

(96) 



,,c,3 . / / ^ n ( ' ' \ (-^"(q^' qa))(-^/3(q;^^ q/3)) ^ , . 

'^2a/3lPaqa;P/3q/3j " (-^FalPa, q^; 2j — ^ 2 7^777^2 ^<-^F/3lP/3, q/?; ^j 

2M^ 2/i^ 

(97) 



and 



It is seen that individual terms in eqs. (94)- (97) can be unambiguously identified with the 
pertinent terms in the multiple-scattering series that emerges from the system of equations 
displayed in section 4.3. Moreover, as already mentioned before, it is seen that, in order to 
make these two series coincide, one should omit the terms of the type 6aG^0p with a ^ (3 in 
the equation for Rq,/3 (third line in eq. (73)). 

The triple scattering is relegated to Appendix E. It can be checked again that to this 
order the multiple-scattering series which are produced by the system of equations displayed 
in section 4.3 can be unambiguously identified. The generalization to the higher orders is 
straight forward . 

4.7 Induced three-body force 

The inclusion of the induced three-body force, which is described by the operator K4, does 
not lead to any complication, since K4 is a smooth function of all momenta. Moreover, at this 
stage it is seen that a genuine three-body force could be included without any further ado. 
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4.8 Analytic continuation and the cusp effect 

In the three-body counterpart of the Liischer formula, the on-shell T-matrix elements below 
threshold are necessarily present, due to the fact that the on-shell projector A defined in eq. (53) 
does not vanish for q^ < 0. Instead, a smooth regulator f^ql/^i^) is introduced in eq. (53), which 
effectively suppresses the contributions with —Qq > fi"^. However, since in the original sum 
the variable p^ was positively defined, a question naturally arises, whether such an analytic 
continuation is needed at all. Indeed, as it is known, one can avoid this procedure in the 
two-particle case, even if multiple scattering channels are considered. 

The reason why the analytic continuation in the three-particle case can not be avoided, lies 
in the following: if, instead of the smooth function /, a sharp cutoff ^(^o) is introduced, the 
principal-value integral will have a unitary cusp at threshold, proportional to the factor y— ^. 
In the two-particle case, the quantity Qq depends on z only and the cusp does not cause a 
problem. In the three-particle case, Qq depends in addition on the spectator momentum p^, see 
eq. (62). Due to the presence of the cusp, one can not apply the regular summation theorem: 
the finite volume corrections are suppressed by a power of L, not by exponentials. In order to 
see this, let us consider a simple example with the cusp present 

^ = 73 E IpI ^(^' - P') = r ^'k |k| -^ Yl ^'(P - k) • (99) 

Using Poisson's summation formula, the above expression can be transformed into 

^ //31, ^_ /-A ^31, 

(^IkK E / (^Ik|e'-^ /../.. (100) 

Carrying out the integration in the second term, it is easy to see that it vanishes as L~^ and 
not as an exponential, see Appendix D. 

To summarize, the price to pay for the exponential fall-off of the finite- volume corrections in 
the regular terms^ that contain principal-value integrals is that these principal-value integrals 
should be defined through the analytic continuation below threshold. If one wishes to use the 
information from the physical region only, the finite-volume corrections in the regular terms 
will be power-suppressed, not exponentially suppressed. 

5 Conclusions 

i) In this paper we have derived the three-body counterpart of the Liischer formula. The 
pertinent expressions are displayed in section 4.3. 

ii) The fundamental property of the finite- volume spectrum, which follows from this formula, 
is that the spectrum is completely determined by the S'-matrix elements for the transitions 



^More precisely, the finite-volume corrections in the regular terms fall off faster than any inverse power of 
L [40]. 
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2 — !■ 2,2 — !■ 3 and 3 — >■ 3 in the infinite volume. Consequently, two different potential 
models with the same S'-matrix elements lead to the same spectra up to the exponentially 
suppressed corrections. 

iii) The equations given in section 4.3 have a complicated structure. However, due to the pres- 
ence of the on-shell factor A(k^, ^q), the dimensionality of equations is reduced as compared 
to the original Faddeev equations. Namely, in the Faddeev equations we have two mo- 
menta describing the three-particle intermediate state. The integration over one of these 
momenta is removed by the on-shell factor. This is similar to the conventional Liischer 
formula, which becomes an algebraic equation (after the truncation of the partial-wave 
expansion), whereas the original Lippmann-Schwinger equation was an integral equation 
in one momentum variable. Despite this simplification, a direct numerical solution of the 
three-body equations in a finite volume, as described in section 2.2, may still prove to be 
less challenging. 

iv) In this paper, we restrict ourselves to the non-relativistic potential model. Considering 
the processes within the field theory will introduce several novel aspects, e.g., relativistic 
effects, particle creation and annihilation, etc. However, we do not expect that these 
effects will upset the proof given in the present paper. Further investigations are planned 
in this direction in the future, and the results will be reported elsewhere. 

v) It is legitimate to ask, whether it is possible to use equations from section 4.3 in order to 
extract the S'-matrix elements from the measured spectrum. Due to the complex nature 
of these equations, a straightforward extraction can be very complicated. Adopting a 
strategy similar to that of ref. [33] - introducing a phenomenologically reasonable param- 
eterization of the potential and fitting the parameters of the potential to the lattice data 
- could be more promising. The physical observables in the infinite volume (S'-matrix ele- 
ments, resonance pole positions) can be then obtained from the solution of the scattering 
equations. 

vi) From the above discussion it becomes clear that it would be very interesting to carry 
out calculations of the finite- volume spectrum in different realistic models (see, e.g. [44, 
45]), which predict both the Roper resonance and iV(1535) in the infinite volume. These 
calculations may shed light on the puzzle related to the "wrong" level ordering for such 
systems in lattice QCD. The findings of the present paper guarantee that the spectrum 
will stay stable with respect to the choice of a model and the variation of its parameters, 
provided the S'-matrix elements in the infinite volume remain the same. 
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A Three identical particles with a separable potential 

In this Appendix, we list the expressions, which appeared first in section 2, in a special case 
of the pair potentials having separable form. Further, we restrict ourselves to the case of three 
identical particles with the masses mi = m2 = rris = m. 

The free-particle states of n identical particles are normalized, according to 

|ki.--k„) = 4^at(ki)---at(k„)|0). (A.l) 



In case of the identical particles, there is only one term in the two-body Hamiltonian H2-S.2, 
given by eq. (1) - the sum over a,/3 drops out. Introducing the appropriate symmetry factors, 
the two-particle potential in the separable model can be written in the form 

v^a(q',q) = ^t'(q')^(q)- (a.2) 

Assuming that the vertex Fq, is also separable 

f.(q';p,q) = ^A(q')/(p,q). (A.3) 

The effective two-particle potential then can be written as 

(q'|w(z)|q) = -t<q>(q) + X{q')d{z)X{q) , (A.4) 

where d{z) is given by 

d^p d^q P{p,q) 



d{z) 



(27r)3 (27r)3 Sp^ q2 

Am m 



dV c?V rf3p d^q /(p',q')(p'q'|m(2;)|pq)/(p,q) 



(2.)= (2.)3 (2.)3 (2.)3 ^3,„ ^ 3(pr ^ (qT _ , _ ,,\ (3„ ^ M + 3! ^ , _ ,0 

4m m / V 4m m 



(A.5) 



where 

3p^ 



(p'qV(^)lpq) = ^ E^'«){(2^)'^'(P^-P/^)^( 



z — m. 



Am. 

■.,13-. ■ 



^ , 3(p^ 

-f T\ z — mi — 



' \2\ / Qti2 



4m 
with 



{p'JY{z)\p,)t(^z - ^ - ^) }^(q/3) , (A.6) 



1 1 

pi = -2P + q' P2 = -2P-q' P3 = p, 

3 1 3 1 ,, ^, 

qi = -^p-2q' q2 = ^p-2q' q3 = q> (A-T) 
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and similarly for p'^,, q^. 

Further, the quantity t{z) is defined by 



^ ^ 2m + — - ;z - iO 



771 

Finally, the quantity (p'|Y(2;)|p) obeys the equation 



(p'|Y(z)|p) = 2Z(p', P) + / ^ 2Z(p', p")r (^z - m - ^^) (p"|Y(^)|p) , (A.9) 



with the kernel 



7(W r.^ ^(pV2 + p)t;(p^ + p/2) 

m 



B Separable potential in a finite volume 

Below, as in Appendix A, we restrict ourselves to the case of three identical particles interacting 
with the separable pair potentials. The finite-volume versions of eqs. (A.5)-(A.9) are obtained 
by merely replacing the integrals by the sums. We do not display these equations explicitly, 
with one exception: the finite-volume version of eq. (A. 8) is given by 



r^(^) = -(l 



q 2m + -^-2;^ 
m 

where, as already mentioned above, the sum over the momentum q runs over (cf. eq. (32)) 

q = ^(l + in), l,n = Z^ (B.2) 

Here, the shift of the discrete variable 1 is related to the CM motion of a two-particle pair in 
the rest frame of three particles. 

Using the regular summation theorem [40], we obtain that, up to the terms exponentially 
suppressed at large values of L, 



77iv'^{p) \ 27r^ Jq q^ — p 



n2 



o-o(i;^-)^;^g ,,,^._,. , ^^t ^--(t 



(B.3) 
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In other words, the components 9i of the vector 6 are either or vr, for the even/odd components 
of the vector n. 

It is interesting to compare the above equations to the three-body equations in a finite 
volume, which were considered in refs. [46-48]. A straightforward comparison leads to the 
conclusion that the sole difference between these equations consists in the replacement Zqq{1\ z/^) 
by Zoo(l; i''^)i i- e., in neglecting the CM motion of the two-particle sub-systems in the rest frame 
of three particles. Note also that the parameter is related to the topological phase considered 
in ref. [49]. 

C Summation in eq. (41) 

The multiple-scattering series in eq. (40) can be reproduced through the following system of 
linear equations 

%\ /1\ /O y2 ysX /TA 
T2 = 1 + h/i 1/3 T2 . (C.l) 

It is straightforward to ensure that 

- (Ti + T2 + T3 - 1) = 1 + (yi + y2 + 2/3) + (2/1(2/2 + 2/3) + 2/2(^1 + 2/3) + 2/3(2/1 + 2/2)) H , 

(C.2) 

which has exactly the same structure as the series in eq. (40). 

On the other hand, eq. (C.l) can be solved directly. The solution of this equation is given 
in eq. (41). 

D The asymptotic behavior of the last term in eq. (100) 

Carrying out the integration in k, the last term in eq. (100) can be written as 

l2 = :r^^ \-T^ + 2ReJ4{x) + 2x lmJs{x)-x''ReJ2ix)\, x = LA , (D.l) 

^ /I xJ I ±1. I 

nGZ\0 ■- I I -> 

where 

neZ\0 ' ' 

Using the relation from ref. [39] 

J2z\^\=gp{z) gp{z) = {es{0,z)f, (D.3) 

nGZ 
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where Osi^O, z) is the elhptic ^-function 

oo 

fc=— 00 
which obeys the following integral representation, 

^3(0, z) = -i / du z"' cot(7ru) , (D.5) 

a set of differential equations, which relate Ji{z) with gp{z), can be derived. For example, for 
i = 1, the pertinent equation has the form 

dJijz) gp{z) - 1 

— - — = . (D.6) 

az z 

Equations for i = 2, 3, 4, ■ ■ ■ can be obtained in the similar fashion. Integrating these differential 
equations, the following expression for Jj(2) can be straightforwardly obtained 

JAz) = T^^ ['dy{lnyy-' ^''^'y^~\ . = 2,3,4. (D.7) 

(« - 1)! Jo y 

Further, from eq. (D.4) it follows that 

\es{0,yen\<Os{0,y), 9s{0,y)>l, fory>0. (D.8) 

From the above equation, we readily obtain 

miO,ye^nf-M<mO,y)f-~l. (D.9) 

Consequently, 

\Mz)\ < j-^ r^i/llnyri^^M^. (D.IO) 

(« - 1)! Jo y 

The above integral converges a.t y = 0. Further, since the series in eq. (D.4) converges for 
\y\ < 1, the divergence in the integral may occur only on the upper limit y = 1. Indeed, using 
the integral representation, one finds that 63(0, y) ~ (1 — y)~^^^ as y — ?■ 1~. However, due to 
the fact that \ny vanishes as y — )■ 1, the singularity in the integrand is of the integrable type. 
Consequently, |Jj(2;)|are uniformly bound from above and therefore, the quantity L^l2 is also 
bound at L — )• 00. 

Last but not least, we wish to address here the issue of using a sharp cutoff at a momentum 
A. As seen, e.g., from eqs. (D.l) and (D.2), the expression for the quantity I2 contains a sum of 
rapidly oscillating terms proportional to sin(|n|LA) and cos(|n|LA), as L — )■ 00. These terms 
are the artifacts of using a sharp cutoff and should disappear when the cutoff is removed, since 
no observable effect in the infrared should emerge from the ultraviolet cutoff. However, the 
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expressions diverge at A — >■ oo. In order to tackle this problem, the simplest way is to introduce 
an additional smooth cutoff in the expressions 

h=Y. ^"(^)' '^"(^)= / 7^ exp(-5p2) |p|e^"P^ (D.ll) 

nez\o "^ ^ '^' 

and then consider the limits A — )■ oo and e — !■ (in this order). It is easy to show that 

1 1 f^ F 

lim lim J„ = — - lim / rfpp^e"^'^ sinp, e' = -^-—r . (D.12) 

e-i-OA^oo 27r2L^ n^ e'^OJo n^L^ 

Now using the equality [50] 

lta^~.ppV=>%i„p= Un,(-A{±,F,(l;|;-±)}) = -2, (D.13) 

where iFi{a; b; z) denotes the conffuent hypergeometric function, we finally get 

In Y. j^. = oiL-') ■ (D-14) 



TT^L^ ^ In 

nGZ\0 ' ' 

Comparing with eq. (D.l), one sees that all the oscillating terms disappear, as expected. 
However, the final expression is still only power-suppressed in L, not exponentially suppressed. 

E Triple scattering diagram 

Below we shall consider the triple scattering diagram, which is given by 



/. 



Z^ ri5 L^ 



f a(q'; vlO{-y.{<. v', + j^^ vD) 



"a 



X 



2M^ 2/i;3 
(f^(q;p^q^))* 



(E.l) 



2M^ 2/x^ y V 2Af^ 2/i^ 
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In this expression, first the summations over the momenta q^, q-^ are carried out, similarly as 
in section 4.6. Further, with the use of the formal relation (see Appendix F) 

(27r)353(p, _ ^^ 



^^ ^ (k'. + ^;^^ V'L? .. .„ .. (K + T7^^ P '' 



2Ma 2/i^ J V 2M^ 2/i^ 



Pa! ^ JG'k7( P7,P/3 



G^Ka ( Pa, P/3 + T7 1 Pa5 ^ 1 G'k7 ( P7, P^ + TF^ P"/^ ^ 



Gf4plp'B + T:r^—p'L-^' 
+ 



""' P.? 



2 ^p;^ + M_^ P7) 



(p^,)^ ^^^ M-m. 
2M7^ 2^^ 



2M/3 2^,^ 2 Ma 2fla 



2M^ 2/i/3 2Mj 2/i^ 



Gf7 I P7, P^ + V7— — P7; ^ 

+ ^ m^ , (E.2) 

M + Ml! + !!!:^ZEK!!1 _ , 

2Ma 2/i„ 

as well as the symmetry properties of the denominators 



(E.3) 
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the above expression can be rewritten in the following form 
2 ^ J (27r)3(27r)3 



Pa 

+ T^^l 1^^^"^^ ' P^<)^3i7(p^<' p;^' P7q7)(f 7(q; P7q7))* 
+ L^l^J (27r)3 (27r)3 "^^ ' Paq^) ^3^/37 (Paq^; P/?; P7q7)(r7(q; P7q7)) 

P7 

+ Z^ E / |^r.(q'; p'^q::)f/3t,,(p'^q::; p;.; P7q7)(f 7(q; P7q7))* 

PJ3P7 

+ ^6 5Z / (^;^r«(q''Paqa)^3i7(Paqa;p^;P7q7)(r7(q;P7q7))* 

+ i^ E /^r.(q'p'^q")^3i7(p'X;p;.;P7q7)(f7(q;P7q7))*}, (e.4) 

PaP7 

where 

uif^M^-^ p^5 P7q7) = gUpI <; ^){-VaK, €))gUpI h'L; z){-Uq';„ q'p) 

X GKfs{p'f},q!p;z){-V^{q!!^,q^))GK^{p^,q.Y,z) (E.5) 

Usip,ipWL; P'fs; P,^,) = G,4pI q'i; z){-VMl q:))GKa{pl q'L; ^)(-^/3(q^, q'^')) 

X GK/sip'fs, q^'; z){-v^{q!^, q7))G'K7(P7' q7; ^) 

+ gUpI q^; z){-v.{ql q':,))GUpl q«; ^)(-^/j(q^, q")) 

X G'K/j(p^,q'/^';2;)(-'l^7(q^,q7))GK7(P7'q7;^) (E-6) 
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gUpI <; ^)i-VMl €))Gkp{p'0, q^; z){-Uq';„ q)^')) 

Gf/jIp^, q$'; 2;)(-V'7(q^, q7))GK7(P7, q^; z) 



X 



X 



+ 



GKaipl q'L; z){-Uq:, q:))G,p{p'^, q;^; z){-Uq% q);')) 



X G'f/3(p^, q)^'; 2;)(-Ky(q", q^))GK7(P7, q^; z] 



(E.7) 



^3;/37(p'X;p^;P7qo 



= GK„(p::,q::;z)(-v;(q'^,C))GK/.(p;3,q;^;^)(-^/3(q;^,q")) 

X Gk7(P7' q'v ^)(-'^7(qr q7))'^F7(P7' q7' ^) 

+ GKaiPa,<^a] z){-Va{ql,q'^))GKp{p'p,^^] z){-Vp{ql^,q^)) 

X G'F7(p7,q^;z)(-'p7(q^,q7))GF7(p7,q7;2:) (E.8) 



Usafi-yiMa^ P^! P7q7) = GKa{Pa, ^^a^ z){-V^{q!^, q'^))G'F/3(p'^, q'^; z) 

(-^/3(q'^q;^'))(-^7(qrq7))^ , , 

■GF7(P7,q7;z) 



M + 



iPn 



2ikL 



+ 



iq: 



//\2 



2/^7 



(E.9) 



t^S7(Paqa;p'/55P7q7) = GF„(p'^,q'I;^) 



-K(q:i,q:;))(-y,(q^,q-)) 



M 



(p;.)^ + M! _ , 



2M;3 2/i^ 
X Gf/3(p^, q^; 2;)(-1/T,(q", q7))G'K7(P7, q7; z 



(E.IO) 
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t^3i,(p«q«; p^; P7q7) = gUpI €; ^){-VaK, €))gUp:, €; z){-v,{^;, q-)) 



*-^FaiPa, qc,; ^j— T / / x2 — /-//n2 ^<-^F7iP7, q^; 2;j (,— K^i^q^, q7jj<j-F7iP7, q7; 2J 



, r< (11 II \i \r I II -ii\\n i " -" x (~^^(q^' ^/^ ))(^^'(q7' q7)) /^ / n 

+ G'PaiPa, q^; 2;)(-K.(q„, q^))GFa(Pa, q^; ^j^ jz^^ T~fJ\2 ^G'f7(P7' ^7' ^) 



^ (p^ lq^_ 
2M^ 2/i^ 



(E.ll) 



where 



-W ^^7 / 

■»" - p^^aF^p- (^■") 

It is straightforward to observe that the above terms reproduce the multiple-scattering series 
of the Faddeev equations in a finite volume (see section 4.3) up to Oiy^. 

F Product of two energy denominators 

Consider the expression 

J = -Pi^ 1—^ 1 T^JV77^5 1 T^ = TTy^0(p)c^ic^2, (F.l) 

L3 ^:-^ (a2 - (p + ci)2)(62 _ (p + c2)2) L^ ^^-^ 

where 0(p) denotes a regular function. 

If Ci = C2 = and a^ 7^ 6^, the proof of the desired relation given in eq. (E.2) immediately 
follows from the identity 

1 ^ { ^ ^\ rF2i 

(a2-p2)(62_p2) 52_a2\^a2_p2 ^2 _ p2 ^ ^'^ 

by using the splitting from eq. (52) in the individual terms on the r.h.s. of this equation. 
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In the generic case, let us consider the partial- wave expansion 

0(p + Ci-Ci) . 0i(p + Ci' 



{p)d2 



62 - (p + Ci - (Ci - C2))2 b'^-{p + Ci- (Ci - C2))2 

Y^yUP + Ci)X,(ci - C2)iriK(p + c^f) , (F.3) 



Irn 

and, analogously, 

0(p)rfl = J2 y^^iP + C2)3^L(Cl - C2)K2li{p + C2)2) . (F.4) 

Im 

Further, define the quantities 

0^ = $^:yim(p + c2):yL(ci-c2)i^2K^') (f.s) 

im 

Then, the quantity J from eq. (F.l) can be rewritten as 
p 



kd2 - Oidi(j)id2 - 92d2((>2di 



+ 9idi(pid2 + 92d2(p2di 



If 



6, = /(a7/x^)^(A^ - (p + cY) , ^2 = f{bVfi')e{A' - (p + C2)^) . (F.6) 

a^ > , 6^ > , 

Ka^ - b^) - (ci - C2)'| > 2\a\ |ci - C2I, |(a' - 6^) + (ci - 02)'! > 2|6| |ci - C2I ,(F.7) 



then the quantities 01^2; 02'^i are non-singular and one may use the regular summation 
theorem. Namely, the expression in the square brackets in eq. (F.6) is non-singular, so the 
summation can be replaced by integration there. Using the same technique as in section 3, 
one straightforwardly arrives at eq. (E.2). Finally, the relation for the generic values of the 
parameters a,b,Ci,C2 can be obtained by analytic continuation of both sides of eq. (E.2) in 
these parameters. 
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